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Several phenomena in phase transition, fracture mechanics, liquid crystals, can be mo- 
delled as energy minimization problems where the natural energy has both volume and 
surface terms. In many cases the energy functional is obtained as a limit of approximating 
functionals and some of its properties can be deduced from the approximation process. 

A basic step is then to obtain, starting from these properties, an integral representation 
of the energy. We consider here this problem for local functionals T defined on the space BD 
of functions with bounded deformation, which are lower semicontinuous with respect to the 
L 1 -topology, satisfy linear growth and coercivity conditions, as set functions are (restrictions 
of) Radon measures, and are invariant with respect to rigid motions. In order to identify 
the volume and the surface densities we follow the global method for relaxation introduced 
by Bouchitte, Fonseca and Mascarenhas in jjj for functionals defined on the space BV of 
functions with bounded variation, which is characterized by the identification of both bulk 
and surface densities from a local Dirichlet problem. This kind of approach has already 
been used in some other contexts, as, for instance, homogenization, where the homogenized 
density is obtained from a Dirichlet problem in the cell. 

An example of functional in the class we consider is given by the relaxed functional F of 
the bulk energy 



with respect to the L 1 -norm topology, where O is a bounded open subset of and / is a 
Borel function satisfying standard linear growth assumptions. The integral representation 
of F on BD(f2) was studied by Barroso, Fonseca and Toader in ||, where the global method 
was applied in order to derive the surface density, while the volume density was obtained 
by a direct proof using the explicit form of the functional F. 

In this paper, the bulk density is deduced from the global method and the approximate 
differentiability of BD functions proved by Ambrosio, Coscia and Dal Maso in [jlj, while the 
surface density is obtained exactly as in || . Note that both our result and the one in || are 
valid for functions in SBD(f2), i.e. integrable functions u for which the Cantor part E c u of 
the measure Eu vanishes. An integral representation in all the space BD(f2) would require 
more information on E c u, since the only property that E c u vanishes on 7i N ~ 1 -a finite Borel 
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subsets, proved in 0], is not sufficient. We recall in Section 2 some useful properties of BD 
functions. 

In Section 3 we prove the integral representation theorem (Theorem [D^) and give an 
example showing why we assume the invariance with respect to rigid motions. In the last 



section we apply Theorem 3.3 to obtain the integral representation in SBD(£1) of some T- 
limits arising in the homogenization of multi-dimensional structures recently studied in the 
context of linear elasticity and perfect plasticity by Ansini and Ebobisse in || , following the 
measure-theoretic approach introduced by Ansini, Braides and Chiado Piat in 0. 

2. Notation and preliminaries 

Let N > 1 be an integer. We denote by M NxN the space of N x N matrices and by 
M-sym the subspace of symmetric matrices in M NxN . For any £ £ M NxN , £ T is the 
transposition of £. Given u, v £ R , u <g> v and «0d := (u ® v + v <E) u)/2 denote the tensor 
and symmetric products of u and v, respectively. We use the standard notation, C N and 
H N ~ 1 to denote respectively the Lebesgue and (N — l)-dimensional Hausdorff measures. 

Let ft be a bounded open subset of 1^; we denote by B(f2), O(0) and Ooo(Q) the family 
of Borel, open and open subsets of f2 with Lipschitz boundary, respectively. For any x £ 
and p > 0, we denote by B(x,p) the open ball of R N centered at x with radius p, by 
Q(x,p) the cube of centre x and sidelength p, while Q v (x,p) is the cube with two its faces 
perpendicular to the unit vector v. When x — and p — 1 we simply write B and Q. If p 
is a Radon measure, we denote by \p\ its total variation. 

Definition 2.1. A function u : fl — > M. N is with bounded deformation in fl ifu £ L^Q, R N ) 
and Eu := (Du + Du T )/2 £ Mb(Q, M^^^), where Du is the distributional gradient of u 
and M(,(n, M^^) is the space of -valued Radon measures with finite total variation 

in f2. 



The space BD(57) of functions with bounded deformation in il, introduced in flijj , has 
been widely studied, for instance by Anzellotti-Giaquinta Kohn pQj , Suquet p2[ , and 
Temam Jl3[ . It is a Banach space when equipped with the norm 

IMIbD(O) : = II m IIl1(J2,R«) + \Eu\(Q). 

It is sometimes convenient to consider also the distance between two functions u, v £ 
BD(ft) given by 

d(u,v) := \\u-v\\ L1{ntRN) + \\Eu\(Q)-\Eu\(Q)\. 
The topology induced by this distance in BD(f2) is called intermediate topology. We denote 

by A the convergence with respect to this topology. It is well known (see Temam |l| ) that 
the trace operator tr : BD(£!) — > L 1 (r, M. N ) is continuous when BD(f2) is equipped with the 
intermediate topology. 

Whenever the open set fl is assumed to be connected, the kernel of the operator E is 
the class of rigid motions denoted here by 1Z, and composed of affine maps of the form 
Mx + b, where M is a skew-symmetric N x TV matrix and b £ R N . Therefore 1Z is closed 
and finite-dimensional. 

Fine properties of BD functions were studied, for instance, in [jD, || and [ fl0| . We recall 
that if u £ BD(f2), then the jump set J u of u is a countably (H N , n — l)-rcctifiable Borel 
set and the following decomposition of the measure Eu holds 

(2.1) Eu = £uC N + E s u = £uC N + ([«] v u )K N ^ 1 \-J u + E c u , 

where [u] := u + — u~ , u + and u~ are the one-sided Lebesgue limits of u with respect to the 
measure theoretic normal v u of J u , £u is the density of the absolutely continuous part of 
Eu with respect to C N , E s u is the singular part, and E c u is the Cantor part and vanishes 
on the Borel sets that are cr-finite with respect to H N ~ 1 (sec [Q). 
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Moreover, the following theorem on the approximate differentiability of BD functions was 
proved in flj . 

Theorem 2.2. Let be a bounded open set in with Lipschitz boundary. Let u G BD(J7). 
Then for C almost every x G fi there exists an N x N matrix Vu(x) such that 

(2.2) lim * [ Hv)-u{x)-Vu{ X ){y-x)\ dy = Q 

P N JQ{x,p) P 

and 

(2.3) lim -1 / \Hy)-u{x)- Eu(x)( y -x),y-x)\ dy = Q 

P JQ(x,p) \v - A 

for C N almost every x G il- 
ia particular, by (|]^) u is approximately differentiable C N almost everywhere in Q and 
the function Vu satisfies the weak L 1 estimate 

t 

where C(N, f2) is a positive constant depending only on N and fi. 
From (2.3) and ( |2.2| ) one can easily see that 

(2.4) £u{x) = (Vu(x) + Vu(i) T )/2 for £ N -&.c. x G Q. 

Analogously to the space SBV introduced by De Giorgi and Ambrosio in || , the space SBD 
was defined in Q. 

Definition 2.3. The space SBD(f2) of special functions with bounded deformation, is the 
space of functions u G BD(f2) such that the measure E c u in j\2.j[) is zero. 



C N ({x G n : |V«(aO| >t}< ||u|| BD(n) Wt > 0, 



3. Main result 

Let T : BD(fi) x Ooo(fl) — > [0, +oo] be a functional satisfying the properties mentioned 
in the introduction, more precisely, 
(f) J-(-,A) is L 1 (A,M. N ) lower semicontinuous; 

(2) for every u G BD(Q), 

(3.1) ^\Eu\(A) < T{ Ul A) < C(C N (A) + \Eu\(A)); 

(3) !F(u, •) is the restriction to Ooo(fl) of a Radon measure; 

(4) T(u + R) — J-(u) for every u G BD(£1) and every rigid motion R. 

Since the properties (2) and (3) give the absolute continuity of T(u, ■) with respect to the 
measure /i := C N + |-E s u|, in order to obtain the integral representation of T , we need only 
to identify the volume and the surface densities whenever u G SBD(f2). To do this we define, 
as in JjJ, see also [||, for every u G BD(f2) and every A G 

m(u, A) := mf{F(v, A) : v G BD(ft) , v\ 9A = u\ 9A }. 

The basic idea of the global method in consists in comparing the asymptotic behaviours 
of m(u,Q(xo,e)) and J-(u, Q(xq, e)) with respect to fi(Q(xo,s)) as e — ► + , and to show 
via a blow-up argument that, the volume and surface densities are obtained from a local 
Dirichlet problem (see Lemma 3.2). 

We shall use the following lemmas, similar to Lemmas 3.1 and 3.5 in [Q for BV-functions, 
proved in the case of BD-functions in H Lemmas 3.10, 3.12]. 

Lemma 3.1. (§, Lemma 3.10]) There exists a positive constant C such that for any U\, 
i*2 G BD(f2) and any A G Ooo(O) we have 

|m(u!,A) -m(u 2 ,A)\ < C [ |tr (it! - u 2 ){x)\dH N - 1 (x). 

JdA 
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Lemma 3.2. (@, Lemma 3.12]) If T satisfies conditions (l)-(3) then 

T{u,Q v (x Ql e)) m(u,Q u (x ,E)) iV _ 1 
lim — — — — = lmi — — -— \i a.e. 6 " ana for all v € b 

e^O fJ,{Q v (X ,£)) e^a fx[Q u {x ,e)) 

We prove now the integral representation result. 

Theorem 3.3. Let T : BD(f2) x Ooo(£l) — » [0, +oo] be a functional satisfying properties 
(l)-(4-). Then for every u £ SBD(f2) and A £ Ooc(fl) we have 

(3.2) F(u,A)= [ f{x,Su)dx + [ g{x, [u^rfdH 1 *- 1 , 

JA J,J(u)nA 

where 

to Q^ -p/ c n r m(^(- - x ),Q(x ,e)) 

(3.3) /(zo,0 := hmsup j, 

(3.4) 9{xo,\, v) := hmsup rr— 

£-+0 £ 

for all i £Sl,Ael N ,(e M^*, i/ G S^" 1 , and w/iere 

u . fy) ._ I A ifyv>Q 

u X ,u{V) — | o otherwise. 
We use the same notation for •) and its extension to the Borel subsets of O. 



Proof, (i) T/ie volume part. Let u G SBD(f2) and choose xo <E £1 such that 
dF(u,-) F(u,Q(x ,e)) 
-dC^ {Xo)= l^o 

(3.6) lim / \ u ( x ) ~ u(xo) — Vu(xo)(x — x^)\dx = 0, 



(3.5) — -j-^ — (xq) = lim - — ^ — ■ exists and is finite, 



(3.7) limp 7 |.Eu|(Q(a:o 1 e)) = |£u(a:o)|, 

(3.8) lim^|E s u|(Q^ 0) £)) = 0, 

(3.9) -5csn(*o) = ^ • 

Let, for every y 6 Q, 

u e (y) := and u (y) := Vu{x )y. 



By (3.6) the functions u £ converge to iio in L^QjR^). Moreover, 

|S« e |(Q) |Buo|(Q). 

Indeed, by definition 

|.Eu e |(Q) = sup / (kv<f>(y)dy 

^(q.m"*") Q 
\\<t>\\oo<l 

= sup — f (u(x) — u(xo))div (p(x)dx 

,€Oj(,(. ,O,C) ^ 
llvl|oo<l 

= T^I^KQOce)), 

where y>(x) := <j){*=p-). Then from Q we get |J5tt e |(Q) -> |£u(x )| = |Ett |(Q), where 
we used also the formula (2.4). This shows that u e — > uq in BD(Q). 
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On the other hand from the continuity of the trace with respect to the intermediate topology 
it follows that 

tr(u e (y)-Vu(so)(tf))|dW* r - 1 (i/) 

dQ 

' 1 \tr(u(x)-u(x )-Vu(x )(x-x ))\dH N - 1 {x) -» 0. 

dQ(x ,e 



Then by (3J3), Lemmas 3J and 3J2 we obtain 

(LF(u,-) m(u,Q(x ,e)) 
__ (a;o) = to _ 

m(«(i ) + Vm(i )(' - s ), Q(io,£)) 

e^O e N 

Now condition (4) with R{x) := u(xq) + Vu ( x o)-^ u ( x o) ^ x _ x ^ i m pij es that 
dT(u,-) m(£u(x Q )(- - x ),Q(x ,e)) 

-dE^ {xa) = H ^ 

= f(x ,£u(x )). 
(ii) The surface part. As in |^|, Proposition 5.1], it can be proved that 



J-(u,APiJ u )= / g(x,u + ,u , 
Jj(u)nA 



v)dU N - 1 

(u)nA 



where 

/ \ a \ v m(u x ,e,v(- - x ),Q u (x ,e)) 
g(x ,\,O,v) := hmsup ^ 

for all x S fi, A, 6 G K w , e S^" 1 , and where 

, > f A if y • v > 
u A ,fl,v(y) := I fl otherwise. 

Using again condition (4) we obtain 

m(uA,0,r/(- - £o), Qr/^o, e)) = m(u A -e^(- - x ), Q u (x ,e)), 

hence ( |3.4[ ), concluding thus the proof. □ 

Remark 3.4. As a particular case the result in || is recovered, i.e. if F is the localized 
lower semicontinuous envelope of the functional F given by (1.1), then 

F(u,A)= [ SQf(£u(x))dx+ [ (SQ^^O^W)^- 1 ^) 

for every u G SBD(O) and every A G O^Cl), where SQf is the symmetric quasiconvex 
envelope of / introduced by Ebobisse in Q , and characterized by 

SQf(0 = inf | £ /(£ + £^(x))dx; <p G W£'°°(A, R^)} , 

for every £ G and for every bounded open subset A of R N , and /°° is the recession 

function of /. 



Remark 3.5. Note that hypothesis (4) is not a consequence of hypotheses (l)-(3). In fact, 
without condition (4) of Theorem [3J^ we would obtain that 

T{u,A)= I f(x,u(x),'S7u(x))dx+ / g(x,u + (x),u~(x),v u (x)) dH N ~ 1 (x) 
J A Jj u nA 
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for every (it, A) e SBD(ft) x Ooo(fi). In particular, for every u G W 1 - 1 ^, R N ), 

J-(u,Q) = I f(x,u(x),Vu(x))dx, 
Jn 

which, under some continuity assumption on T with respect to u, for instance, assuming 
that there exists a modulus of continuity ip(t) satisfying 

\T(u(- —z) + w,z + A) — T{u, A)\ < ip(\w\ + \z\)(C N (A) + \Eu\{A)) , 

for all (u,A,w,z) E BD(fi) x Coo(fi) x R N x R N , such that z + AcO, implies that for 
C N -almost every xo £ ^ and every p G R N , the function f(xo,p, •) is quasiconvex. By (2), 

<f(x ,tx ,t) <c(i+\z+e\)- 

The following example shows that there exists a rank-one convex function <f) '■ M 2x2 — ► 
[0, +oo[ which satisfies 

(3.io) ^ie+e T i <0(O <c(i + ie + e T i) veeM 2x2 , 

and which depends also on the antisymmetric part of the matrix £. Let £ := ^ ^ ^ ^ . It 

is enough to define such a function on the matrix f ^ ^ ^ and then to add the quantity 

i ui c- /0 6\ c + 6/0 1 \ c-6/0 -1 \ . ., 
|a| + |a|. Since I ^ / = — 2 — I 1 / 2 — 1 1 j ' W6 & * unc * 10n 

/i(6, c) which is separately convex, satisfies the linear growth condition 

< h(b,c) < C(l + \b + c\) 
and which depends also on c — b. An example of such a function is the following: 

— (x + y) if x + y < 

if0<a; + y<l and xy < 

xy if x, y > and max(a;, y) < 1 

x + y — 1 if a; + y > 1 and max(x, y) > 1. 
Therefore, the function </> : M 2x2 — » [0, +00 [ given by 

a b 



h(x,y) := < 



</>(£) : = h(b, c) + |a| + \d\ + \b + c| where £ = 



c 



is rank-one convex, has the linear growth (3.1C), and does not depend only on the symmetric 
part of the matrix £. 

4. Application: homogenization of periodic multi-dimensional structures 
In jj| , the authors studied the asymptotic behaviour of functionals of the form 

F e (u,Q) := J <p(^~, -^-) d ^e 

defined on a particular class of functions with bounded deformation, denoted by LDP c (Q), 
and given by the functions u € L p (£l,R N ) whose deformation tensor Eu is an absolutely 
continuous measure with respect to n £ , with p-summable density dEu/dfi e , where n £ is 
defined by ^ £ (B) := e N n{\B), with fi a fixed 1-periodic Radon measure and tp is a Borel 
function 1-periodic in the first variable. Assuming the standard p-growth condition on if and 
that the measure /1 is 'p-homogenizable' (see |], §4]), the authors proved a homogenization 
theorem (Theorem 5.1). Precisely, they proved the existence of the T-limit Fh om of the 
functionals F e with respect to L p -convergence in the Sobolev space W 1,p (Sl, R N ), and with 
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respect to L 1 -convergence in BD(f2). They showed that the T-limit admits the integral 
representation 

(4.1) F hom (u, O) = / tp h om{Eu)dx 
in W 1,p (f2; M. N ); moreover, if ip is convex and p = 1 then 

F h 

in BD(f2), where tphom is described by an asymptotic formula. However, in the second 
case, the question about the integral representation of the T-limit without the convexity 
assumption on tp remained open. Notice that such an assumption is too strong in the 
vectorial calculus of variations. As shown in ||, (see the proof of Theorem 5.1), the T-limit 
verifies the properties (l)-(3) and the invariance with respect to rigid motions follows from 
the fact that the approximating functionals F e have this property. So we can apply Theorem 
|3.3| to obtain that 

(4.2) F hom (u,A)= [ f(x,£u)dx+ f g(x, [u], u^dU*' 1 , 

J A Jj u nA 

for every u G SBD(Sl) and every A G Ooo(fl). Now, from the integral representation (4.1) 
in W 1 ' 1 (Q,'§l N ) and the relaxation theorem 3.5 in [El, one can easily see that 

f(x,0 = ¥W(£) and g(x, a, v) = p™ om {a © v) 

for every x G il, a G M. N , v G S N ~ X , and for every £ G M^^. Notice that, since </3hom is 
symmetric quasiconvex, that is 



<£hom(£) < f V'hom^ + £l/j{x))dx 
J A 

for every tp G W 1,oo (A, R^), £ G and for every bounded open subset A of R N , then 

^hom 



Acknowledgements. The authors wish to thank G. Dal Maso for many useful dis- 
cussions concerning the subject of this paper, in particular for suggesting the example in 



Remark 3.5. The work of Rodica Toader is part of the European Research Training Network 



"Homogenization and Multiple Scales" under contract HPRN-2000-00109. 

References 



[i 

[2 
[3 
[i 
[5 
[6 

[7: 

[8 
[9 
[10 



Ambrosio L., A. Coscia and G. Dal Maso: Fine properties of functions with bounded deformation. Arch. 
Rat. Mech. Anal. 139 (1997), 201-238. 

Ansini N., A. Braides and V. Chiado Piat: Homogenization of periodic multi-dimensional structures. 
Boll. Un. Mat. Ital. (8) 2-B (1999), 735-758. 

Ansini N. and F. B. Ebobisse: Homogenization of periodic multi-dimensional structures II: the linearly 
elastic/perfect plastic case. Adv. Math. Set. Appl. To appear. 

Anzellotti G. and M. Giaquinta: Existence of the displacement field for an elasto-plastic body subject 
to Henky's law and Von Mises' yield condition. Manuscripta Math. 32 (1980), 101-131. 
Barroso A. C, I. Fonseca and R. Toader: A relaxation theorem in the space of functions of bounded 
deformation. Ann. Sc. Norm. Sup. Pisa 29 (2000), 19-49. 

Bellettini G., A. Coscia and G. Dal Maso: Special functions of bounded deformation. Math. Z. 228 
(1998), 337-351. 

Bouchitte G., I. Fonseca and L. Mascarenhas: A global method for relaxation. Arch. Rat. Mech. Anal. 
145 (1998), 45-68. 

De Giorgi E. and L. Ambrosio: Un nuovo tipo di funzionale del calcolo delle variazioni. Atti Accad. 
Naz. Lincei Rend. CI. Sci. Mat. Natur. 82 (1988), 199-210. 

Ebobisse F.: On lower semicontinuity of integral functionals in LD(Q). Ric. di Matematica XLIX 
(2000), 65-76. 

Kohn R. V.: New Estimates for Deformations in terms of their Strains, Ph.D. Thesis, Princeton Uni- 
versity, 1979. 



8 



FRANCOIS EBOBISSE AND RODICA TOADER 



[11] Matthies H., G. Strang and E. Christiansen: The Saddle Point of a Differential Program. Energy 

Methods in Finite Element Analysis, Wiley, New York, 1979. 
[12] Suquet P.M.: Un espace fonctionnel pour les equations de la plasticite. Ann. Fac. Sci. Toulouse 1 

(1979), 77-87. 

[13] Temam R. : Problemes Mathematiques en Plasticite. Gauthier-Villars, 1983. 

(Frangois Ebobisse) S.I.S.S.A., Via Beirut 2-4, 34014, Trieste, Italy 
E-mail address, F. Ebobisse: ebobisse@sissa.it 



(Rodica Toader) S.I.S.S.A., Via Beirut 2-4, 34014, Trieste, Italy 
E-mail address, R. Toader: toaderOsissa.it 



